Emergence of Particle-Hole Symmetry near Optimal Doping in High- Temperature 

Copper Oxide Superconductors 
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High-temperature copper oxide superconductors (cuprates) display unconventional physics when 
they are lightly doped whereas the standard theory of metals prevails in the opposite regime. For 
example, the thermoelectric power, that is the voltage that develops across a sample in response to 
a temperature gradient, changes sign abruptly near optimal doping in a wide class of cuprates, a 
stark departure from the standard theory of metals in which the thermopower vanishes only when 
one electron exists per site. We show that this effect arises from proximity to a state in which 
particle-hole symmetry is dynamically generated. The operative mechanism is dynamical spectral 
weight transfer from states that lie at least 2eV away from the chemical potential. The emergence of 
this symmetry close to optimal doping points to pairing in the cuprates being driven by high-energy 
electronic states. 



I. INTRODUCTION 



As charge carriers are doped into the wide family 
of parent copper-oxide ceramics (cuprates), supercon- 
ductivity obtains with a maximum transition tempera- 
ture at a particular chemical composition. The result 
is a dome-shaped superconducting region in the doping- 
temperature plane. The origin of the dome remains piv- 
otal to the solution of this problem because the prop- 
erties of all high T c superconductors change drastically 
around optimal doping, x c . Unconventional physics per- 
vades for x < x c in which the standard theory of metals 
breaks down whereas such a description is recovered in 
the opposite regime. The efficient cause underlying this 
drastic change in the physics of the cuprates as the dome 
is traversed remains one of the key mysteries of these ma- 
terials. To solve this problem, it is instructive to focus 
on a correlate of superconductivity, that is a phenomenon 
which supervenes on superconductivity and exhibits an 
abrupt signature at the top of the dome. Such a corre- 
late of superconductivity could elucidate the ubiquitous 
dome-like structure in the doping-temperature plane of 
the cuprates. 

To this end, we consider the thermo-electric power (S), 
defined as the voltage difference that develops across a 
material in response to a temperature gradient. Phys- 
ically, the thermopower measures the entropy per car- 
rier and is a direct probe of particle number rather than 

I 



charge conservation, a distinction which is of utmost 
importance in this work. A universal feature observed 
amongst different families of hole-doped cuprates^— is 
that the thermo-electric power is non-zero everywhere 
except at one particular doping value. Shown in Fig. 
(QJi) is a collation of thermopower data at 290K for a 
wide range of hole-doped cuprates as a function of the 
planar hole density per unit area, P p i, which measures 
the hole concentration in each Cu02 layer. P p \ pro- 
vides an even-handed way of comparing different families 
of cuprates that is independent of the doping method 
or sample quality, be they single crystals or not. As 
long as x refers to the number of holes per Cu02 layer, 
x and P p i can be used interchangeably as we do here. 
Clearly shown is that for all hole-doped cuprates, the 
thermopower at 290K is positive for P p \ < 0.24, whereas 
for P p i > 0.24, it is negative. To correlate this sign 
change with superconductivity, Honma and Hor— used 
the thermopower doping scale as opposed to the Pres- 
land scale 6 , 1 - P c /T c max = 82.6(z - 0.16) 2 , in which the 
maximum T c is artificially fixed to be 0.16, to determine 
the doping level at which T c is maximized. We show in 
Fig. §\)p) the maximum superconducting transition tem- 
perature T™ ax as a function of P p i for 23 different types 
of cuprates classified by the number of copper-oxide lay- 
ers in each unit cell. Of the 23 cuprates shown, only four 
single-layer materials exhibit any significant deviation of 
the maximum T c from the planar hole density at which 
the thermopower vanishes. 



One might argue that this agreement is irrelevant since the thermopower is highly temperature dependent. However, 
experiments 7 indicate that at 110K, the zero-crossing of the thermopower is shifted to 0.26, quite similar to the value 
at 290K . In addition, one might also question the validity of the thermopower scale to calibrate the doping level. To 
address this question, we reprint here a figure (see Fig. ([5])) from their paper which illustrates that the thermopower 
scale is in excellent agreement with the doping level in Y123 determined by three different experimental methods. 
Hence, the zero-crossing in the vicinity of T™ ax is a robust feature of the cuprates that has direct bearing on why T c is 
so high. Our analysis indicates that mixing with states that lie at least 2eV (the optical gap in the cuprates ^ 1:L ' 13 i 14 ) 
away from the chemical potential acount for the dramatic sign change of the thermopower and hence are the cause of 
the high transition temperature for most of the cuprates, particularly those exhibiting T c > 70K (see Fig. (QJl). 
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FIG. 1. a) Universal behavior of the thermoelectric powei>i 
(290K) as a function of planar hole density (P p i), for var- 
ious families of hole-doped cuprates. All exhibit a sign 
change at P p \ — 0.23. Above the solid-bold horizontal 
line, the thermpower obeys the functional form, S 290 (P p i) = 
392exp(-19.7P p i) for 0.01 < P pl < 0.21. Below the solid-bold 
horizontal line, S 290 (P p i) = 40.47 - 163.4P p i for 0.21 < P p[ < 
0.34. These functional forms were used^ to determine the 
hole doping levels for all the cuprates rather than the widely 
used empirical formula 6 1 - T c /T™^ = 82.6(:r - 0.16) 2 which 
artificially fixes the optimal doping level of all cuprates to be 
0.16. b) Maximum transition temperature as a function of the 
planar hole density using the thermopower scale to determine 
the doping level. Except for three single-layer materials, the 
vanishing of the thermopower coincides with the doping level 
at which the transition temperature is maximized. 



II. THERMOPOWER AND SPECTRAL 
WEIGHT TRANSFER 



The dramatic sign change of the thermopower in the 
cuprates poses a distinct theoretical problem because it 
represents a stark departure from the predictions of free- 
electron physics. For example, in the standard theory of 
metals, Fermi liquid theory, the sign of the thermopower 
unambiguously reveals the sign of the dominant charge 
carriers in the material. Consequently, in the standard 
theory of metals, the thermopower only vanishes in a 
half-filled band where the number of filled and empty 
states is equal, the particle-hole symmetric condition. In 
the cuprates, the parent materials possess a half-filled 
band. Doping away from half-filling by the introduction 
of holes should only increase the particle-hole asymme- 
try. Consequently, it is unexpected that the thermopower 



should vanish at x = 0.24 corresponding to a filling of 
n = 0.76. Existing treatments of the thermopower in 
the cuprates address either the low-doping regime^ or 
invoke special band structure effects such as Van Hove 
singularities^ and the shape of the Fermi surface^ to 
describe the sign change. However, given that the sign 
change obtains near optimal doping, the approximate 
crossover from strong to weakly interacting physics, the 
explanation of the sign change of S, must incorporate the 
physics of strong correlations. 

To address this question, we focus on the minimal 
model that captures the strong correlations of the copper- 
oxide planes of the cuprates, namely the one-band Hub- 
bard model, 
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FIG. 2. Hole-doping level from various techniques compared with the doping scale extracted from the thermopower, P p i. The 
red points are obtained from (near-edge X-ray absorption fine structure) NEXAFS: red circles are OD-Y123 1 ^, red diamonds 
are co-doped Ca-YaC123 (at x=0.1}^, red squares are Calcium-doped Y1236^. The blue points are from nuclear-quadrupole 
resonance (NQR) measurements: left blue arrow is OD-Y123^- and the right blue arrow is OD-Ti2201 1 . The green points are 
from ARPES: green circles are Strontium-doped La214^, green squares are also Strontium-doped La214^, uprrow, overdoped 
BiPb2201 19 , down arrow co-doped BiP2201— . The star corresponds to angular magnetoresistance oscillations (AMRO) 20 . 
Reprinted from Phys. Rev. B 77, 84520 (2008). 



are nearest neighbours, Ci„ annihilates an electron with 
spin a on lattice site i, t is the nearest- neighbour hop- 
ping matrix element and U the energy cost when two 
electrons doubly occupy the same site. Our conclusions 
carry over naturally to any n-band model of the cuprates 
as long as the largest energy scale is the on-site en- 
ergy, U in Eq. (QJ. For the cuprates U ~ A.OeV and 
t ~ O.beV. Consequently, the cuprates reside in the 
strong-coupling regime where standard weak perturba- 
tive treatments breakdown. 

Quite generally, the strong correlations in the Hub- 
bard model mediate a universal sign change in the ther- 
mopower at non-integer fillings as seen in the cuprates. 
Since U 3> t in the cuprates, perturbation around the 
atomic limit (vanishing hopping, t = 0) rather than the 
non- interacting regime (U = 0) is the correct starting 
point. In 1974, Beni3^ considered this limit of the Hub- 
bard model and calculated the exact expression for the 
thermopower 
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to 0(t/U) and U 3> fcgT, where kg is Boltzmann's con- 
stant and e is the electric charge. This expression van- 



ishes exactly at a hole doping level of x = 1/3. However, 
the reason why it vanishes at x = 1/3 has never been 
understood. Consider the general expressio n 25 ' 26 for the 
thermopower, 
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with the transport integrals Lij in the relaxation-time 
approximation given by 



L; 






and t{ui) being the relaxation time, 
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In the atomic limit, the single-particle spectral function, 
A(k,u>), which when summed over momentum yields the 
density of states, has no momentum dependence. Con- 
sequently, from the form of L12, a vanishing of the ther- 
mopower arises entirely particle-hole symmetry. How 
does such a particle- hole symmetry obtain for x = 1/3 
in a Hubbard system? Alternatively, what is so special 



about the (2a;/ 1 — x) ratio in the argument of the log- 
arithm in Eq. ([5])? The answer lies in spectral weight 
transfer which was only understood as the signature 
physical phenomeno n 27 i 29 i 31 of strongly correlated sys- 
tems in the early 1990's, roughly 20 years after Eq. ([2]) 
was derived. As no one yet has presented an explana- 
tion of the Beni2i result in this light, we offer one here. 
In the atomic limit of the Hubbard model, as depicted 
in Fig. (J3j) , one particle resides per site. In a Hubbard 
system consisting of N sites, there are N ways to re- 
move an electron and N ways to add one, constituting 
the photoemission and inverse photoemission bands, re- 
spectively. In the atomic limit, the splitting between 
these bands (also known as the lower and upper Hub- 
bard bands) is the Mott gap, U, as shown clearly in Fig. 
Q. Now consider removing a particle. There are now 
N — 1 ways to remove a particle and N + 1 ways to add 
a particle. However, only N — 1 such states lie in the 
high energy scale as there are only N — 1 ways to add 
a particle to the system so that it costs an energy U. 
The two remaining states correspond to the two ways of 
adding an electron with either spin up or spin down to 
the empty site. Such addition processes cost no energy 
and hence must lie directly above the chemical poten- 
tial. In general if x holes are introduced, the spectral 
weight immediately above the chemical potential grows 
as 2x whereas the weight in the lower band decreases as 
1 — x. The ratio in the Beni21 expression is simply the 
number of states above and below the chemical potential 
in the photoemission band in the atomic limit. When 
the two spectral weights are equal, particle-hole symme- 
try obtains and the thermopower vanishes. The exact 
particle-hole symmetry condition for the photoemission 
band of a Mott insulator is 2x = 1 — x. As a result, in a 
strongly-correlated system, spectral weight transfer and 
the thermopower are intimately linked. 

Real Mott systems are not in the atomic limit because 
t/U ^ 0. This results in two new effects, the second 
of which destroys the necessity for a strict particle-hole 
symmetry to be present for the thermopower to vanish. 
First, the number of particle addition states immediately 
above the chemical potential, 
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L= I N(uj)duj, 



(6) 



computed by integrating the single-particle density of 
states from the chemical potential to a cut-off energy that 
demarcates the low from the high energy scale, exceeds 
2a?2&. The reason is simple. Away from the atomic limit, 
no eigenstate of the Hubbard model involves a fixed num- 
ber of doubly occupied sites. Harris and Lange^ that 
such mixing with the high-energy scale leads to an in- 
crease in the weight of the lower band to 1 + x + a, where 
a represents all the correction due to the mixing with the 
upper band, a is strictly positive as shown by Harris and 
Lange^ and confirmed by numerics 2 ^. Because the total 
weight of both bands is 2, the weight of the upper band is 
l — x — a. In general a is doping dependent and can be as 
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FIG. 3. Evolution of the single-particle density of states from 
half-filling to the one-hole limit in a doped Mott insulator 
described by the Hubbard model. Removal of an electron 
results in two empty states at low energy as opposed to one 
in the band- insulator limit. The key difference with the Fermi 
liquid is that the total weight spectral weight carried by the 
lower Hubbard band (analogue of the valence band in a Fermi 
liquid) is not a constant but a function of the filling. 



large as 20%2£. To interpret the thermopower in terms 
of the entropy, it is helpful to reinstate an atomic-limit 
picture of the LHB. We have showr>2£ how this can be 
done by re-defining doping level to include the holes that 
are dynamically generated as a result of mixing with the 
upper Hubbard band. To this end, we define x' = x + a. 
Consequently, the filled part of the spectrum at low ener- 
gies has a weight of 1 — a;' and the empty part 2a;'. Hence, 
L = 2x' > 2a; as a result of dynamical spectral weight 
transfer. The new particle-hole symmetric condition for 
the lower band is 1 — x — a = 2(x + a) or a; pris = 1/3 — a 
which is strictly less than 1/3, the atomic-limit value. 
However, a strict particle-hole symmetry is not neces- 
sary for the thermopower to vanish once the hopping is 
turned on. At play here is the fact that the spectral func- 
tion is momentum dependent for t ^ 0. As a result, the 
vanishing of Lyi no longer arises from a simple balancing 
between the density of states above and below the chem- 
ical potential, but rather t{lo) averaged above and below 
the chemical potential must be equal. This pushes the 
zero crossing of the thermopower to even lower doping 
levels. We note importantly that computational tech- 
niques will differ in the spectral function which in turn 
will affect the magnitude of the thermopower. However, 
the existence of the zero crossing and its location relative 
to x < 1/3 — a is a robust feature determined entirely 
by the integrated weights of the spectral intensity below 
and above the chemical potential and not on the compu- 
tational method used. 



III. REALISTIC CALCULATIONS 

To quantify the arguments presented here, we com- 
puted Eq. ((3} for the Hubbard model using the cellular- 
dynamical mean- field theory 32 . In this method, we con- 
sider exactly the dynamics on a 4-site (plaquette) cluster 
and treat the rest of the lattice as a bath whose degrees 
of freedom are integrated over. The coupling of the clus- 
ter to the bath is thus treated in a mean-field fashion and 
the cluster quantities are evaluated self consistently using 
the cumulant lattice reconstruction scheme 3 ^. While no 
variational principle exists for this method, several tests 
in limiting cases 3 ^ in which the answers are known ex- 
actly have revealed that the local correlations that the 
CDMFT method builds in provide an adequate starting 
point for quantifying the physics of strong correlations. 

Using the CDMFT method, we obtained the one- 
particle spectral function A(k, ui) for various values of 
hole doping x and on-site Coulomb repulsion U (ex- 
pressed in units of the nearest neighbor hopping integral 
t) . From Fig. (g]), we see that both U = 8t and U = At 
at T = O.li (slightly higher than the 290K of the ex- 
periments) exhibit a sign change before the atomic limit 
value of x — 1/3. For U — At, the sign change occurs at 
x s« 0.19 whereas for U — 8t, S vanishes at x = 0.21 in 
agreement with the perturbative argument that as t/U 
decreases, the critical value of the doping at which the 
thermopower changes sign must increase. For U = 8t, 
both the magnitude of the thermopower and the S = 
condition are in agreement with the experimental values 
in Fig. (QJi). The inset shows the difference between the 
spectral weight above and below the chemical potential 
for U — 8t. As is evident, this quantity is finite even 
when the thermopower vanishes, in direct constrast to a 
non-interacting system in which particle-hole symmetry 
is essential for S = 0. The deviation from the particle- 
hole symmetric point, which lies in close proximity to the 
doping level at which S vanishes, is due entirely to the 
momentum dependence of the spectral function at finite 
t/U . Consequently, the vanishing of the thermopower 
in the cuprates at a doping level significantly less than 
the atomic limit of x — 1/3 is a signature that the dy- 
namical contributions (the t/U part) to the low-energy 
spectral weight and the momentum dependence of the 
spectral function cannot be ignored. Such an extreme 
sensitivity of the thermopower to the spectral weight re- 
distribution obtains because the thermopower is gener- 
ated by a thermal gradient, unlike an electrical gradient 
in a Hall measurement. Consequently, the thermopower 
is fundamentally an experiment about the conservation 
of particle number (as opposed to charge conservation 
in the Hall experiment), which is of course governed by 
electron spectral weight shifts. 

Based on the calculations presented here and our anal- 
ysis of Eq. ([2]), we correlate the vanishing of the ther- 
mopower at the top of the dome for a wide range of 
cuprates to a proximity to the particle-hole symmetric 
state which is driven predominantly by mixing (dynam- 
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FIG. 4. Thermoelectric power computed for two different 
values of U=8t, plotted as a function of hole doping, x. The 
solid line is Eq. ([2]), the exact result for the thermopower of a 
doped Mott insulator in the atomic limit. In the inset A is the 
difference of the spectral weight above and below the chemical 
potential. As is evident, S vanishes at a doping level slightly 
lower than that at which A = 0. In a strongly correlated 
system, naive particle-hole symmetry is not essential for the 
vanishing of the thermopower. 



ical spectral weight transfer) with states in the upper 
Hubbard band. As it is the dynamical part of L that 
is relevant here, maximising T c is tied to states at high 
energy, in particular states that lie 2eV away from the 
chemical potential. This conclusion is consistent with the 
results from optical experiment a 34 ' 35 which indicate that 
the condensation energy in the superconducting state 
arises from electronic states that lie 2eV away from the 
chemical potential. 

Aside from influencing the vanishing of the ther- 
mopower, dynamical spectral weight transfer also pro- 
vides a general mechanism 36 for the breakdown of Fermi 
liquid theory in doped Mott insulators in the absence of 
any symmetry breaking. In a Fermi liquid, the phase 
space available to add a particle at low energies (that is 
L) is exhausted by enumerating the number of ways (rih) 
of adding an electron to the empty states. This principle 
follows from the Fermi liquid tenet that at any chemical 
potential, the number of quasiparticle excitations equals 
the number of bare electrons in the system. This ac- 
counting clearly breaks down in a real doped Mott sys- 
tem because the number of ways of adding an electron is 
simply 2x but the phase space available to add a particle 
at low energies, namely L, exceeds^ 2x as shown above. 
Hence, the number of low-energy particle addition states 
per electron per spin exceeds unity. As a result in a doped 
Mott insulator, there are electronic states at low energy 
that have no counterpart in the non-interacting system 
(L/rih > 1) and Fermi liquid theory breaks down. 



IV. FINAL REMARKS 



We propose here that close to the top of the 
dome where particle-hole symmetry is dynamically 
generate d 37 ! 38 , a quantum critical point obtains sig- 
nalling a transition to weak-coupling physics. That is, 
in the strongly overdoped regime, the coupling to the 
high-energy scale ceases, giving rise to L/rih = 1, and 
traditional descriptions in terms of Landau quasiparticles 
apply. Recent soft x-ray Oxygen K-edge experiments 39 
indicate that L/nu does saturate to a doping indepen- 
dent value in the overdoped regime once the pseudogap 
terminates as predicted here. Similar experiments should 
be performed as a function of temperature below and 
above the T* line. According to the theory we have re- 



cently constructed 3 ^ to explain the onset of non-Fermi 
liquid behaviour in a doped Mott insulator, L/n^ be- 
low T* should increase. The increase should be given 
by a, for which quantitative estimates have been ob- 
tained recently 3 -. The fundamental problem of super- 
conductivity in the cuprates reduces to uncloaking pre- 
cisely how the high-energy scale physics mediated by dy- 
namical spectral weight transfer creates a phase coherent 
condensate. 
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